Monitoring D-region Variability from Lightning Measurements 
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Abstract 

In situ measurements of ionospheric D-region characteristics are somewhat scarce and rely mostly on sounding 
rockets. Remote sensing techniques employing Very Low Frequency (VLF) transmitters can provide electron density 
estimates from subionospheric wave propagation modeling. Here we discuss how lightning waveform measurements, 
namely sferics and tweeks, can be used for monitoring the D-region variability and day-night transition, and for local 
electron density estimates. A brief comparison among D-region aeronomy models is also presented. 


1. Introduction 

The ionospheric D-region is extremely difficult to study not only because the chemistry is complex but in situ 
measurements are scarce and sometimes unreliable. Unlike the regions above where only positive ions and electrons are 
important, negative ions also play a key role in the D-region. Additionally, presence of metallic particles, provided by 
meteoritic bombardment, and heavy cluster ions produce intricate chemical reactions in the mesosphere. On the other 
hand, the D-region altitude range is too low for satellites and too high for atmospheric balloons. Except for sounding 
rocket in situ measurements made mostly during day, the D-region is usually investigated through propagation 
techniques that effectively provide data at much larger scales, employing satellites, rockets, and ground-based 
transmitters. An alternative to Very Low Frequency (VLF) ground-based transmitters is using lightning signatures, 
namely sferics and tweeks, to study the D-region variability through subionospheric wave propagation modeling. In this 
work we report VLF electric field measurements made onboard a stratospheric balloon during the African Monsoon 
Multidisciplinary Analyses (AMMA) SCOUT campaign [1, 2]. The balloon flight lasted about 3 h, encompassing the 
time of sunset, thus contributing to assessing the variation of the propagation conditions due to the D-region electron 
density changing. After a brief description of the experiment and atmospheric conditions of the flight, we provide the 
results of an analysis of several hundred impulsive electromagnetic events for which the transverse resonance mode of 
the surface-ionosphere waveguide is positively identified. These observations are in fair agreement with a wave 
propagation model that considers an electron density distribution near dusk derived from the International Reference 
Ionosphere (IRI) model [3], considering three representations of the chemical processes. 


2. Instrument Description and Balloon Measurements 

The experiment took place on 7 August 2006 in Niamey, Niger (geographic coordinates 13. 5N, 2E), and the 
balloon was launched at TO =1645 UT. With an ascent velocity of -5 ms' 1 , the balloon reached the ceiling altitude of 23 
km at 1802 UT, traveling -200 km almost exactly westward until the end of the flight. During the ceiling period, the 
balloon altitude decreased slowly down to about 20 km after 1845 UT; the last measurements were recorded at 19 km. 
The sunset at the balloon position occurred at 1843 UT (TB, 118 min from TO); the sunset at the position of the 
strongest local convective system (-400 km east of the balloon final position) and in the D-region above occurred 93 



(TS) and 130 min (TD) from TO, respectively. Images from the AMMA campaign meteorological radar showed active, 
though variable, convective systems about -100 km east and south of Niamey [2], Numerous lightning and strong 
electric field transients were observed. 

The HV-AIRS (High Voltage - Atmospheric Impact of Radiations and Sprites) electric field experiment was part 
of a multiple instrument flight chain that included other sensors devoted to observations of aerosols and atmospheric 
composition. It measures the atmospheric electric field from DC to -4 kHz employing the double probe technique, 
which has been widely used onboard balloons and satellites to measure electric fields. The effective length of the 
antenna is about 40 cm and the DC and AC components of the electric field are measured in two independent channels. 
The AC signals are measured in the frequency range from 4 Hz to 4 kHz, the noise level is 30 mVm' 1 Hz' 1/2 , and 
multiple gain settings allow for measuring electric fields up to ±2.5 Vm" 1 . Although aiming a proof-of-concept for 
planetary missions - consequently, the physical parameters are very much limited - instrument sensitivity and sampling 
are sufficient for sferics and tweeks detection. The HV-AIRS gondola also includes simple optical detectors to identify 
lightning as far as -300 km after sunset and a GPS receiver. Figure 1 illustrates typical sferics and tweek waveforms 
recorded during the flight. Most waveforms are detected during nighttime because propagation conditions are more 
favorable, i.e., the ionosphere is sharper and offers better reflection conditions for VLF subionospheric wave 
propagation. 
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Figure 1: Typical (left) day and (middle) night sferics recorded by HV-AIRS onboard the stratospheric balloon at 
altitude of about 20 km. (right) Distinctive tweek recorded during nighttime at 19 km (adapted from [6]). 


3. Numerical Modeling 

To understand the balloon data, we have developed a wave propagation finite element model that benefits from 
algorithms previously used to study wave propagation on Earth and other planetary environments [4, 5]. The model is 
run in the transient and harmonic wave propagation modes to compute the electric fields in the surface-ionosphere 
waveguide. For the transient mode, we consider a vertical current element with a typical lightning current profile. The 
stationary mode computes the electric field spectrum of the waveguide in the range 0-10 kHz, considering a vertical 
Hertz dipole. Below 60 km, a scalar conductivity profile is used. Above 60 km, we use the parallel, Pedersen, and Hall 
components of the conductivity tensor, considering the Appleton-Hartree equation which describes the refractive index 
in cold magnetized plasma [6]. The electron density is derived from the International Reference Ionosphere (IRI) 
model, including the three parameterizations of the D-region: (i) default (IRI-2007) - IRI- 1 [7], (ii) Friedrich (FIRI) - 
1RI-2 [8, 9], and Danilov - IRI-3 [10]. Figure 2 shows transient signals at the balloon ceiling altitude computed with the 
transient method. The discharge starts at t = 1 ms and the electric field is computed in several locations, i.e., for various 
angular separations between the radiating source and the balloon. The frequency of the waveform is -2 kHz and the 
propagation delay increases with angular separation. It is worth mentioning that the quasi-static component decreases 
with distance and a clean sinusoidal waveform is observed. Although lightning stroke accurate characterization 
considering intensity, duration, distance, and direction requires triangulation and receiver synchronization, it is possible 
to estimate the source distance from waveform morphology, i.e., shape, attenuation, and quasi-static component 
contribution. Comparison between the balloon data and model waveforms suggests that about 75% of the recorded 
lightning occurred no farther than 350 km, i.e., 3° of angular separation, which is significant for the day-night transition 
accurate detection. 
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Figure 2: Waveforms of a lightning discharge calculated with the transient mode algorithm for several angular 
separations from the source. 


4. Results and Discussion 

From the few thousand impulsive events recorded by the AC channel we are able to retrieve frequency 
information from about 700 sferics and a few tweeks. Three approaches are tested to compute the waveform mean 
period, namely the method of least squares applied to damped sinusoidal functions, Fourier analysis, and crest-crest 
period counting. From typical events shown in Figure 1, we can see that the frequency,/, and wave attenuation, a, are 
lower for nighttime propagation. These results are consistent with wave reflection in a sharper layer at higher altitude, 
in line with IRI predictions. Figure 3 presents the frequencies derived from sferics and tweeks waveforms. Although the 
rate of lightning strokes is roughly constant during flight, the events from which frequency calculation is reliable are 
concentrated after sunset, an indication of more favorable propagation conditions during nighttime. The histogram on 
the right-hand-side provides the number of events as a function of frequency, showing a broad peak about 2 kHz. 



Figure 3: Comparison between frequencies derived from sferics and tweeks and our numerical model, (left) Squares 
correspond to the frequency of individual impulsive events; the solid line and error bars represent the mean value and 
standard deviation calculated over periods of 30 min. The colored lines are numerical results using IRI electron density 
distributions, considering (blue) IRI- 1 , (green) IRI-2, and (red) IRI-3 model options for the IRI D-region. The labels T s , 
T b , and T d represent the sunset time at the strongest convective cell, balloon position, and the D-region above the 
convective cell, respectively, (right) Histogram of impulsive events as function of the mean frequency measured during 
the flight. 




During the solar eclipse of 7 March 1970, Reeve and Rycroft [11] characterize tweeks resulting from nighttime 
propagation conditions to determine the height of the ionospheric reflector, and are able to infer the change of the 
effective height of the ionosphere during the eclipse. They obtain 69 and 76 km at the start/fmish and maximum totality 
of the eclipse. A previous study of the present data shows that the day-night ionospheric depletion occurring in the 
altitude range from 73 to 82 km before and after sunset, respectively [6]. These altitudes are somewhat lower that those 
suggested by other studies [cf. 6 and references therein] because the present data are recorded close to sunset, when the 
ionosphere has not yet reached an ionization minimum. Several conclusions can be drawn from Figure 3: (i) 1RI-2 and 
IRI-3 show sharper day-night transitions than that of IR1-1, in line with the balloon measurements; (ii) the frequency 
range derived from the IRI-3 model shows a better agreement with data despite a 30 min offset; (iii) Although 
consistently underestimating the average frequency by -0.15 kHz, the IRI-2 model reasonably matches the profile. 
Under the assumption these results are representative and can be generalized to other latitudes and longitudes, a more 
elaborate model would perhaps combine IRI-2 electron density magnitude with IRI-3 electron density gradients. 

In this work we show that lightning is a valuable tool to study the day-night transition of the ionospheric D- 
region. A straightforward comparison among three models of the D-region (IRI-1, IRI-2, and IRI-3) unveils their 
relative accuracy for deriving electron density distributions. Careful monitoring of the ionospheric day-night variation 
may eventually provide additional information for investigating the chemical processes taking place in the D-region. 
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